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ABSTRACT: This paper concerns the queuing system M x /G 1, B /l, to which the customers are assumed 
to arrive in batches of random size X according to a compound Poisson process and also are served in 
batches. As soon as the system becomes empty, the server leaves for a vacation of random length V. If 
no customers are available for service after returning from that vacation, the server keeps on taking 
vacations till he finds at least one customer in the queue, then immediately begins to serve the 
customers up to the service capacity B. If more than B customers are present when the server returns 
from a vacation, the first B customers are taken into service. If fewer than B customers are present, all 
waiting customers go into service. Late arrivals are not allowed to join the ongoing service. The steady 
state behavior of this queuing system is derived by an analytic approach to study the queue size 
distribution at a random point as well as a departure point of time under multiple vacation policy. It 
may be noted that the results in [5] and [7] can be obtained as special cases from the results in this 
paper. 
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I. Introduction 

Most queuing models assume that customers are arrived singly and are served singly. But this 
assumption is far from truth when we consider those numerous real-world situations in which customers are 
arrived and are served in batches. We call such queues “batch arrival and batch service queues” (or bulk 
queues). Such queues frequently occur, for example, in loading and unloading of cargoes at a sea-port, 

in traffic signal systems, in a mass transportation system, in digital communication, in computer 
network and in production/ inventory systems (e.g., Takagi [1991] and Doshi [1986, 1990a]). 

Server vacation models are useful for the systems in which the server wants to utilize the idle time for 
different purposes. The major general result for vacation models is the stochastic decomposition result, which 
allows the system to be analyzed by considering separately the distribution of the system size with no vacations 
and the additional queue size due to vacations. This important result was first established by Fuhrmann and 
Cooper [1985] for generalized vacation as well as multiple vacation models, where the server keeps on taking a 
sequence of vacations of random length till it finds at least one unit in the system to start each busy period for 
the M/G/l queuing systems. Later Doshi [1986] extended this result for the single vacation model through a 
different approach, where the server takes exactly one vacation at the end of each busy period. In this model if 
the server finds no units after returning from a vacation, he stays in the system waiting for a unit to arrive. 
Shanthikumar [1988] showed that the queue size decomposition holds even for the M/G/l models with bulk 
arrival, reneging, balking etc. In terms of unfinished work in the system, Boxma and Groenendijk [1987] proved 
the decomposition result for the M/G/l type vacation models. Recently Doshi [1990b] and Leung [1992] 
extended the results of Boxma and Groenendijk [1987], Harris and Marchal [1988] and Shathikumar [1988], 
proved the stochastic decomposition result for unfinished work in the system and additional delay due to the 
vacation times respectively in more general setting. 

At present, however, most studies are devoted to batch arrival queues with vacation because of its 
interdisciplinary character. Considerable efforts have devoted to study these models by Baba [1986], Lee and 
Srinivasan [1989], Lee et.all [1994,1995], Borthakur and Choudhury [1997] of M x /G/1 type queuing models of 
this nature have been served by Chae and Lee [1995] and Medhi [1997]. 

In this paper, we first study the steady state behavior of the queue size distribution for the M X /G 1,B /1 
queue with multiple vacation policy at the stationary point of time as well as departure point of time through an 
analytical approach. Also, we show that the departure point queue size distribution of this model can be 
expressed as convolution of the distributions of three independent random variables, one of which is the queue 
size of the standard M X /G 1,B /1 queue without vacations. Efforts have also been made to interpret other two 
random variables properly. 
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II. M x /G 1,b /1 Queuing System With Multiple Vacation 

In this chapter we consider a bulk queuing system with following specifications: 

1. As soon as a service is finished, the server takes into service a maximum of B customers, (B is called the 
service capacity). 

2. If fewer than B customers are present at a service completion epoch, all the existing customers are taken into 
service. 

3. Customers who arrive while a group of customers are being served cannot join the service even if there is 
space available. 

4. The server takes a vacation of random length V if there are nosremotsuc to serve. 

5. If no customers are available for service after the server returns from a vacation, he keeps on taking vacations 
till he finds at least one customer in the queue and begins to serve. 

6. If the server finds at least one customer waiting for service when he returns from the vacation, the server 
immediately begins to serve the customers up to the capacity B. 

We first write the system state equations for its stationary (random) queue size (including those in service, 
if any) distribution by treating elapsed service time and vacation time as supplementary variables. We now 
define the following notations and probabilities: 

A = arrival rate of batch , 

= service rate , 

B = service capacity , 

P = traffic intensity , 

X = group size random variable of arrival , 

C(£=p( x=k ) ; k>i , 

X(Z) = PGF of X , 

S = service time random variable , 

V = vacation time random variable , 



S(x),V(x) = the probability distribution functions of “S”, “V” , 

s(x),v(X) = the probability density functions of “S”, “V” , 

S x (s ) p V x (s) = laplace stieltjes transforms of “S(x)’\ “V(x)” , 

S° = elapsed service time for the customer in service , 

= elapsed vacation time for the server on a vacation , 

N Q (t) = queue size at time t 

Further, it is assumed that F(0) = 0, V( oo) — 1,5(0) — 0 and 
s (oo) — 1 and that V (x) and S (x) are continuous at x=0, so that 

v{x ) dx=dV{x)f[ 1 - V(x)] 

and 

six') dx=dSix)/[ 1 - 500] 

are the first order differential functions of “V” and “S” respectively. 

Let x Q (t) be the queue size at time “t” and s u (t) be the elapsed service time at time “t”. Further, we 

assume that V°(t) is the elapsed vacation time at time “t”. Let us now introduce the following random 
variable 




0, if the server is on vacation at time "t" 

1, if the server is busy at time " t" 
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So that the supplementary variables 5° CO and F°(t) are introduced in order to obtain a bivariate Markov 



process [NQ(t) f L(t)}, where L(t) = V°(t) ,f Y(t) = 0 and L(t) = s 0 (£)ify(t) = i. 
We define 

PamGO dx = jyrob\N q(X) = = V°(t),x < 

V 0 (t) < x + dx] ; x > 0,m > 0 

p ln (x)dx = lim t ^prob{N Q {t) = n,L(t) = S°{t) f x < 

s°(t) ^ x + rfx] ; x > 0^ ii ^ 1 



Now the analysis of this queuing process at the stationary point of time can be done by using forward 
kolmogorov equations, which under the steady state conditions can be written as: 

Gta) Po + U + *M]po,oM = 0 ; x > 0 , (id 

(£?) Po ^ ^ + U + v (x)]p 0iU U) = i c fe p 0j n-fe W ; 

x > 0, n > 1 , ( 2 . 2 ) 

(£;) p ^»W + U + s(x)]p ln (x) = ; 

x > 0, n > 1 , ( 2 . 3 ) 

Po,o(0) = J 0 ” vUlPo.oM rf* + Sf =1 /" sWpii(ar) dx . 



(2.4) 



Where 



Po,o = J 0 “po,oW d* 

These equations are to be solved under the following boundary conditions, at x=0: 

(2.5) 



Po,o(0) = 


- 2p 0j0 . 








Po,n(0) - 


= 0 ,n > 1 


9 






Pl,n(0) = 


: JT rWpo,„W rfz + J 0 “ 


S'WPj 


and the normalizing condition 








yco f” 

Zj?i=0 jq 


p 0 „(x)dx + 2 


loo 

J Tl= 1 , 


Jo p i^ 


l (jc)dx 


Let us define the 


following PGFS 








Po (*; z) 


= Z"=o z "Po,n 


(x) 


; M 


< 1, 


Po(0 ;z) 


= SS > =0 z "P0,n 


(0) 


; |z| 


< 1, 


Pi(*;z) 


= 2“=lZ"Pi,„ 


M 


; 1*1 


< 1, 


Pi (0; z) 


= Z“ =1 z" PljJl 


(0) 


; 1*1 


< 1 



( 2 . 6 ) 



(2.7) 



( 2 . 8 ) 



0 jj 

Now multiplying equations (2.1) by Z and (2.2) by Z and summing over possible values of n yield. 

Po^W = Hn=o -U + v(*)]^“p 0jW (*) + 

^ H^=i Ek=i tfc^Po^-fe G) , 

From PGFS, we get 
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7- Po (*, z) = -[A + t?(x)]P 0 U, z ) + 

ax 

AZ”=iZfe=iZ n C fe p 0j n-feW , 

Now 

Z?=iZfc=iZ n C k p 0 jll -fc(x) = ZCiP 0j0 W + 

zLi* 2 c fePo J z-feW + 

ZLl^ 3 CfcP0 J 3-fcW + *" 

= 2CiP 0j0 W + 2 2 C 1 p 0j i W + Z 2 C z p 0|0 W + 2 3 C 1 p 0j 2 U) 

+ Z 3 C2Po j iW+Z 3 C 3 Po j oW + ■" 

= ZC 1 {p 0 j 0 (x) + zp 0 j i(x) + Z Z p 0j2 W + Z 3 Po 3 (x) + •■■ } 

+ z 2 c 2 {p 0 j 0 W+zp 0jl W + —}+•*• 
=Zfe=i2 k C fc Zn=0 Z "P0 J nW = Zfe=i2 fe C fe P 0 (x,z). 



Let 



*U) Zfe=i2 fe c fe , 



then 



Z^iZLl^^fePo.n-feW *(Z)P 0 (X,Z), 



then 



dx 



P 0 (x, z) = -[A + v(x)]P 0 (x, z) + AX (z)P 0 (x, z) 



= [WX) - A - t>(x)]P 0 (x,z) , 
dP 0 (x,z) = [ Ax ( z ) _ x _ v (x)]dx 

Integrating with respect to x 

/nP 0 ( x, z) = AX(z)x — Ax — J * v(t)dt + constant 

P 0 (x, z) = A. e xx ^ x . e~ Ax . exp{— J* t?(t) dt} . 

To find A put x=0 

P 0 (O,z) = A = £~ =0 z"p 0 ,„(0) 

= Po,o(°) + 2n=iZ"Po, n ( 0 ) 

Then 

A Poo(0) Ap 0>0 . 

Then we have 

P 0 (x, z) = Ap 00 exp{- f* v(t ) dt}. e“*“ I(z ) ]l ,x > 0 , 



we have that 



v(x) dx = 



dV(:r) 
[ 1 -^ 00 ] ’ 



s(x) dx — 



(2.9) 

dS (x) 
[1-S(x>] 



I IJMER I ISSN: 2249-6645 I 



www.ijmer.com 



I Vol. 5 I Iss.2l Feb. 2015 I 321 



Analysis of a Batch Arrival and Batch Service Queuing System with Multiple Vacations 



Then 

v(x)[l - v(x)] - v ( x ) _ v (x)V(x ) = 

ax ax 

d +v(x) V (x) = v (x) . 

This is a linear differential equation of the first order and to solve it multiply both sides by € " 

^ v ( x )e } o v(t ' )dt V(x) = v(x)e^ v( - t)dt 



— \v(x)e S » v(tUt ] -- 

dx\- J 

by integrating both sides, we get 

v(x)e% H,Ut = J* 



-- v(x)e ! « v{, > dt , 
v(u)e^ l,< ,> dt du 



— gjfl ^(0 ^ ^ 



Then 

v(x) = i - e~ s <> v(t > dt 



Substituting this into equation (2.9), we have 

P 0 Uz) = Ap 0 , 0 [l-KW]e- J[1 -* )l1 ,* > 0 

Then 

Po (z) = j; Po (x,z)dx 

= J” Ap 0 o [1 - dx 

= Ap 0j0 [1 - v(x)]e~ A[1 ~ x{z)]x dx , 

integrating by parts, we get 

, , , l-V*\A-AX(z)\ 

P 0 (z) = Apo.o ' 

P (y, = „ i-y*[A-^w] 

PoUJ Po,0 [i_x(z)] 

Similarly, multiplying equation (2.3) by Z and summing over all values of n we obtain 

j- P a ( x , z) = -[A + s(x>] Pi (x, z) 

+ A S” = i Z J =1 c k p ljn _ k (x)z 71 
= -[A + s(x)] P a (x, z) + AX(z ) P a (x,z) 
= [AX(z) - A - s(x)] P a (x, z) 



jfv(t) rff 
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d Pi (_x,z) _ _ A _ 5 (x)] dx , 

Pi( X,z) 

by integrating with respect to x 

lnP 1 ( x , z ) = AX(z)x — Ax — J * s(t ) dt + constant 

Pi ( x , z ) = B. e AXi - z ^ x . e~ Ax . exp{— s(i) dt] 

= B. e -A[l-X(z)]x eX p £_ j x ^ 

= B[1 - S(x)].e~ A[1 - x ^ ]x 

To find B put x=0 

Pi(0 ,z) = B , 

then 

P i(x,z) = Pi (0, z) [1 - S(x)]. e-m-x^x , x > 0 

(2.10) 

Now 

Pi (0, z) - £” =1 z n p ln (0) 

From equation (2.7) 

= S"=l Z U Iff v(x)P0,n (X) dx 

+ / 0 ”*WPl J n+B (*)<**] 

= ff[ln=lZ n P0,n(x)]v(x)dx 
+ /(,“ [E"= 1 Pl.n+B (X)]S (x)dx 

= /“[P 0 (X z) - p 0 0 (x)]v(x)dx 

+ -jg / 0 °°E”= 1 z1l+B Pi >n +B(x)]s(x)dx 
= /o°°[ P o^ z ) - P 00 ( x )]v( x )dx 

+ ~ *B / 0 " [Pi ( x > x) - Sf=i z* p u (x)]s (x)dx 
= P 0 (x, z)v(x)dx - J" p 00 (x)v(x)dx 
+ r P 1 ( x * z)s(x)dx 

- ^Z?=i zi C p uW 5 

From equation (2.4) 

= J 0 “ PoU z)v{x)dx 

— [2p 0j o - 2?=i J 0 °° PijUM*)***] 

+ ^ r P 1 ( x > z ) s ( x )dx 
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- £?= i zi C Pi^ to s(x)dx 
= Jq Po (x,z)v(x)dx - Ap 0>0 

+ ^ c p i ^ x * z ) s ( x ) dx 

+ ^I,?=i(z B ~ z l ) p u (x) s(x)dx 

From equations (2.9) and (2.10) 

- J" A Poo [l - 

+ ^ /" P a (0,z) [1 5(x)]e _jl[1 “ x(z)la; s(x) dx 

— 2p 0j0 

+ ^£?=i(z* - 

+ ^ C Pi (0, *) e- A[1 - x(2)]x dS(x) 

— Ap 0>0 

+ ^Zf=i(z B ~ z l ) /“ p u (x) s(x)dx 
= Ap 0j0 V*[A - AX(z)\ 

+ Pj. (0, z)S* [A — AX (z)] 

- ^Po.o + ^l B =i(z B - z l ) ff p u (x) s(x)dx 

p 

Multiply by Z to get 

Pitt) ,z)S*[A - AX (z)] - z^ttU) 

= z B Ap 00 - z B Ap Q 0 V*[A - AX(z)] 

- 2f=! (z B - z f ) J“ Pi j( (x) s(x)dx 



Then 



Pi(0,z) 



.4jz B p 0 0 {l-F“ [Jt-AZ (z'jV-Yf^Cz 8 -z l ) p lti (x)s(x)dx 

S'[A-AX(z)]-z B 



Therefore, we get 

Pi to = J” Pi to z) dx 

Substituting from equation (2.10) to obtain 

Pi W = C Pi(0,z) [1 - SWle-*-®' dx 
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_ Az g p 0 0 (1 - y * [A- AX (z) ]} f =1 (z B -z l ) J D ° : p lit jx)s(x)dx 

S*[A-AX(z)]-z b 

J 0 °°[l — 5(*:)] dx 

Integrating by parts we get 



_ pz B Po 0 {1— ^(z)]}— S^ 1 (g g — z 1 ) \"f 



Pi (z) = [ 



S*[A-1X(z)]-z b 



1 



r i-s*U-JixW]-, 

L A-AX(z) * 



pz^Po Qd-r^U-jjfCz)]}] r i-5 x U-jJC(z)], 
_ I- S*U-ZX(zV-z B J L A-AX(z) J 



J*[JL-iUr(z)]-z 
-ff /• B J-i f™ 



_ r £^i(z g -g I )fo PtiU)j(a-)ria; l l-S*U-AA(z)] , 
[ 5 x U-1A'U)]-z f J ^ A-AX(z) ^ 



Then 



p ( \ = z a PM{l-V*[A-AX(zy]}{l-S*[A-AX(.z)]} 
1 W ~~ {l-Z(z)}{S*U-AZ(z)]-z s } 



(l-S*U-ilAr(z)]}2f =1 (z fl -z l ) f 0 °° pi^OO-sCx) dx 
{A-AX(z)}{S*U-AX(z)]-z B } 

Using the normalizing condition (2.8) and taking the limit of [P 0 (z)+P 1 (z)] asZ 



By L’Hopital s rule, we get 



PoW = l 

p o(z) = ^WV-V«M] 



-X\(z) 



at z = 1 



P 0 (l) = -Ar'( 0 )p „, 0 = A£(P)Po , 0 



(i) 



Also, 



Pid) = 



By L’Hopital s rule, we get 



P !<X)U — ~ By applying L’Hopital s rule another time, we get 



at z = 1 



P ri -) _ 2i 2 p QO r'(Q)^'.‘(i)s-yo) 

1 2 (l) [Ax\ (l)S * Ho) + B} 

2^V(1)5^(0) f 0 °° puCaj jQQrf* 

2A AT\(l){ AX\(l)S* \ (0)+B) 

_ ^ 2 po j0 g(F)g(Jf)F(^ 

~~ {-AE(X)E(S)+B} 



, we get 
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then 

Pi CD 



E(S~) i) Jq 00 Pl,i (x) s (x) dx 

A 2 p 0j0 ^aO-PCJO£’W+^C^) Sf = i(B-i) fo' Pi i(x)s(x)dx 

B —AE(X)E(S'} 



Now using the normalizing condition to obtain Pq q as follows: 

P 0 (D + P 1 ( 1 )= l 



Afi(lOp 0j0 



>i 2 p 0 ' 0 E(yyEUOE<Ls'>+E(s-y xf-jfg-O Pj. ^(xys(xydx 



B — A E (_X~) E (S') 

= 1 

A i:t\s ) Po o [ li — A e ( x ) i: ( .v ) ] i a- r> [>t> e 0-0 e ( x ) /- ( -V ) i 

EC^j:f = ^CEt — O JT I>x,i 0*0 sCx'ydx = B — AEtX^ECS') 
ABEiV^o o + BOS') Y>i=W B — O sCx)dx = 

B — AEiX')E{S') 



Then 

Po.o 

Let 

P = 

Then 



_ B-AEOOE(S')-E(S) £f =1 (B-i) /“ V± £x)s(x)dx 
~ ABEiV) 

AE{X)E(.sy 



_ B (1- p)- E (S) sf= i(S - £) / 0 P a ^x^s^xydx 
Po '° _ ABEtV) 

Let P 0 (z) + (x) = P (x) be the PGF of the stationary queue size distribution of this M x /G'' 

B /1 multiple vacation model, then 

pm = poo -r ™ 1 



z B PO,Q(l-V*tt-AXQr)]}{l-S*[A-AX(z)]} 

{1 -X{z)}{S* [A -AX(_z)]-z B } 

{1-5* [A-AXCz)]} Zf=i(z B -z 1 } jf pijOdsix) dx 
[A - AX(z)} [S*[A-AX(z)]-z b } 



+ 



{b (1 - />)- E(S) Sf= x (B - i) / 0 “ Pl dx}{l-V * [A—AX(.z)]} 

ABE(Y) [1 -X(z)] 
ri z G {±-S*[A-AX(zn± 

L± + {S~[A-AXCzy]-z B } J 



{ 1 - 5 *u -AX O)] } Zf = i (Z B -Z l ) J 0 °° pi,i(x)sW dx 

{A-AX(z)}{S*[A-AX(z)']-z B } 



Then 

P(z) = 



[b (l-p) -£(5) Sf =1 (B- 0 PiiWsWd^l-ZJS^ {A-AX{zy\ 

S x [A-AX(z)]-z B 



'-BE{ynA-AX(_z}y 
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_ S| LiO g ~-g ) Sp p l i (x)s(x)dx ^ i 

L A-AX(z) S*[A-AX(z)]-z E J ' 

SPECIAL CASE:(I) 

When B=l, then 

_ (l-pXl-z)S* [A-AXtzVil-V* [A-AXjzm 
{S * [A -A X {»] -z) E(_ V) [A -A X (z)] 
p(z) = r (X-p)(.X-z)S" [A—AXQzy] ±—y* [A-AxCzj]^ 

{ S*[A-AX(_zy]-z } ^ E(YXA-AXCz)]^ 

= P(M x /G/l;z)c(z) 

Where 

P(M X /G/ l;z) — The PGF of the stationary queue size 
distribution of the standard M x /G/1 

queue without vacation. This is well known Pollaczek-Khinchine formula for 
M x /G/± queue. 

— r Ci -/>) Ci - z)s l [a-a x Q)] -| 

{S * [A AX(z)] -z} -* ’ 

and $(Z) — The PGF of the number of units that arrive during the residual life of the vacation 

time 

— r l-V* [A-AXiz)]^ 

E(V)iA-AX(z)V ’ 



Differentiating equation (2.2.1 1) w.r.t. Z and taking the limit as Z —• ► 1, we get 
L = ciP(z) \ 

S k=l 

By applying L’Hopital s rule four of times, we get 

, = a 2 Poo e(_x)e(v 2 ) 





f AEOO[A 2 E 2 {X)E($ 2 ')+A E(S ) e(x Z —Xj— p B<B — 1)] 

2B(l-p) 2 " + 



2ABpE(yj+A 2 pE(XjE(v 2 ) 

2(1— p) 



} 



~(_x)s(_x)dx 

+ 2B 2 (1— p) 2 

E {S') Sf =1 (B (R- 1)- i(i- 1)) Sq p l i {x)s{x)dx ^ 1 ^ 

, _ • \L.YL) 

SPECIAL CASE: (II) 



Substituting with value Po.O into equation (2.12), then put B=l, we get 




aeCxSe iv R s + p + 



a 2 e 2 (x)e(s 2 )+ae(s'je(x 2 

2(1 — p) 




Where E ( V R ") = E ( V 2 ) /2E(F) is the mean residual vacation time and L$ is the mean of the 
stationary queue size distribution of this model. 

III. Departure Point Queue Size Distribution 
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